Distribution of exchange energy in a bond-alternating 5 = 1 quantum spin chain 
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The quasi-one-dimensional bond-alternating S = 1 quantum antiferromagnet [Ni(A r , A r '-bis(3- 
aminopropyl)propane-l,3-diamine(/i-N02)]C104 (NTENP) is studied by single crystal inelastic neu- 
tron scattering. Parameters of the measured dispersion relation for magnetic excitations are com- 
pared to existing numerical results and used to determine the magnitude of bond-strength alterna- 
tion. The measured neutron scattering intensities are also analyzed using the lst-moment sum rules 
for the magnetic dynamic structure factor, to directly determine the modulation of ground state 
exchange energies. These independently determined modulation parameters characterize the level 
of spin dimerization in NTENP. First- principle DMRG calculations are used to study the relation 
between these two quantities. 

PACS numbers: 75.10.Pq,75.40.Gb,75.40.Mg,75.30.Ds 



I. INTRODUCTION 

Integral antiferromagnetic (AF) spin chains are best 
known for having an exotic spin liquid ground state 
with a characteristic gap in the magnetic excitation 
spectrumi^ The Haldane gap has been subject to inten- 
sive theoretical and experimental studies for the last two 
decades, and is by now very well characterized and un- 
derstood. The spin wave function of the Haldane ground 
state is not known exactly, but is similar to the easy 
to visualize Valence Bond Solid (VBS) stated The lat- 
ter is constructed by representing each 5 = 1 spin as 
two separate S = 1/2 spins, binding pairs of these into 
antiferromagnetic dimers for each exchange bond, and 
projecting the resulting state back onto the subspace 
where 5 2 = 2 on each site. This wave function is the 
exact ground state of the Affieck-Kennedy-Lieb-Tasaki 
(AKLT) models and is schematically shown in the left 
inset of Fig. ^ Each exchange link carries exactly one 
valence bond, and the periodicity of the underlying crys- 
tal lattice remains intact. Considerably less attention 
has been given to a different quantum spin liquid ground 
state that is realized in integral spin chains with alter- 
nating exchange interactions. As the alternation param- 
eter 5 = (Ji — J2)/(Ji + J-i) deviates from zero (uniform 
chain), the energy gap A decreases and closes at some 
critical value \b\ = S c rj 0.26 ! ^ ( Pi7i9i ! i ) , 1 9i 11 as illustrated 
in Fig. ^ Further increasing \8\ beyond this quantum- 
critical point re-opens the spin gap. The ground state is 
then no longer the Haldane state, but, instead, a dimcr- 
ized one. The corresponding valence bond wave function 
is shown in right inset in Fig. ^ It contains two valence 
bonds on each strong link and none at all on the weaker 
ones. 

The two gapped quantum phases differ by their "hid- 
den" symmetries AiS The Haldane state is characterized 



by a non-vanishing expectation value for the antiferro- 
magnetic string order parameter^ related to a breaking 
of a non-local Z% x Z2 symmetry^ This order parame- 
ter vanishes in the dimerized phase, where Z^ x Zi re- 
mains completely intact. AA 2 - The highly non-local multi- 
spin correlation function that defines antiferromagnetic 
strings can not be expressed through the usual pair spin 



correlation functions (5,- (0)(5j a ''(f)). As a result, the 
"hidden" string order can not be directly observed in 
scattering or other type of experiments. In fact, one does 
not expect any qualitative differences in observable spin 
correlation functions of alternating 5 = 1 chains with 
similar gap energies on different sides of the phase di- 
agram. Distinguishing the two phases in a real 5=1 
alternating-chain compound is therefore a challenging 
task, involving a careful quantitative analysis of the data. 

A model alternating 5 = 1-chain material suitable for 
experimental studies was discovered only recently^ This 
compound is [Ni(A, A'-bis(3-aminopropyl)propane-l,3- 
diamine( / Lt-N0 2 )]C10 4 , NTENP for short, is structurally 
similar to well-known Haldane-gap systems NENP— and 
NDMAPiLii Unlike the latter, NTENP features a dis- 
tinctive alternation of short and long bonds in the anti- 
ferromagnetic 5 = 1 Ni 2+ chains,— Extensive bulk mea- 
surements on NTENP were performed by Narumi et al. 
and are reported in Ref. ^] Susceptibility datai^i^ shows 
that in NTENP A/ J RJ 0.4, almost the same as in a Hal- 
dane spin chainflS Anisotropy effects aside, on the simpli- 
fied phase diagram of an alternating 5 = 1 chainX shown 
in Fig. ^ NTENP must be located somewhere on the 
dashed horizontal line. This line crosses the theoretical 
curve for A (S) twice: at 5 < 5 C (near the Haldane point 
5 = 0) and at S rj 0.40 > 5 C . Bulk measurements can not 
directly probe the microscopic alternation parameter 5. 
Nevertheless, Hagiwara et al. were able to conclude that 
NTENP is in the dimerized phase based on indirect ev- 
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FIG. 1: Field dependence of energy gap in a bond-alternating 
5=1 quantum antiferromaget, based on numerical results by 
S. Yamamoto, Ref. 0- A quantum critical point (dashed line) 
separates the Haldane phase, similar to the VBS state (left 
inset) from the dimerized state (right inset). The possible 
location of NTENP is indicated by the thin horizontal line. 



idence, namely the behavior of non-magnetic impurities 
in this material^ The main purpose of the present work 
is use a microscopic probe (inelastic neutron scattering) 
to directly measure exchange alternation and other cru- 
cial parameters of undoped NTENP. Our experimental 
findings are discussed in comparison with first-principles 
DMRG calculations. 



II. EXPERIMENTAL APPROACH 

Before describing the actual neutron scattering exper- 
iments performed as part of this study, we shall discuss 
the measurement strategies. In particular, we need to 
identify those measurable physical quantities that are 
most sensitive to the effect of bond alternation. 



Structural consideration 




FIG. 2: A schematic view of the bond-alternating S = 1 Ni 2+ 
chains in the triclinic crystal structure of NTENP. The equa- 
torial vertices of the Ni 2+ coordination octahedra are nitrogen 
atoms. The octahedra are coupled via chemically disordered 
NO2 groups. Oxygen and nitrogen atoms are shown as dark 
grey spheres. Light grey spheres are carbon atoms. Hydro- 
gen atoms and intercalated CIO4 solvent molecules are not 
shown. The alternation of short and long bonds in the chains 
is characterized by bond vectors di and d2 , respectively. 



The triclinic crystal structure of NTENP is visualized 
in Fig. E l 14l2 fl The S = 1 chains are composed of Ni 2+ 
ions octahedrally coordinated in an organic environment. 
The chains run along the a axis, one chain per unit cell. 
While all Ni 2+ sites are crystallographically equivalent, 
the Ni-Ni distances within the chains alternate between 
di = 4.28 A and d 2 = 4.86 A. Intra-chain Ni-Ni links 
are covalent and pass trough structurally disordered NO2 
groups. Inter-chain interactions are of Van der Vaals na- 
ture and therefore much weaker. The crystallographic 
symmetry is low, space group PI. The lattice constants 
at room temperature are a = 10.75 A, b — 9.41 A, 



In the following discussion we will use the following co- 
ordinate system: x is chosen along the a* axis, z is along 
c, and y completes a right-handed set of axes. 



B. Model Hamiltonian and observable energy 
scales 

A model Hamiltionian in for NTENP was proposed in 
Ref. In the present work we shall employ a slightly 
different, but equivalent formtiS 
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H = J(l -S)J2 S 2j S 2j+ i + J(l + S)J2 M«-i + } ) 2 



(1) 



r 



Inter-chain spin interactions are expected to be very weak 
and are not included in the above expression. The ex- 
change constant J is, itself, not directly observable ex- 
perimentally. However, it can be reliably inferred from 
an analysis of the measured temperature dependence of 
bulk susceptibility: for NTENP J w 3.4 meVM^ 

A directly measurable quantity related to J is zone- 
boundary energy hu> z h, defined as the minimal energy of 
magnetic excitations with momentum transfer g z b = — , 
where a is the distance between next-nearest-neighbor 
spins. In NTENP a is simply the lattice constant. For 
a uniform (Haldane) chain with 6 = numerical simula- 
tions show that TiUzb ~ 2.7J. 2122 ' 23,2 — It is easy to verify 
that Tiujzb = 2 J for isolated dimers (5 = 1). The com- 
plete 5-dependence of hui z b was recently determined in 
a systematic numerical study^ It is important to note 
that unlike A/ J, Hlo^/J is a monotonic function of 5, 
and can be used to unambiguously determine whether a 
particular material in in the dimerized or Haldane phases. 

Another directly observable energy scale is spin wave 
velocity v. Numerically, v ks 2.49J for a uniform chairi^i 
and v = for the other limiting case of isolated 5 = 1 
dimers. The gap energies A Q for different spin polariza- 
tions are also experimentally accessible. For an isotropic 
uniform chain A w 0.41Ji 21 i 22 i 23 i 24 For isolated dimers 
with 8 = 1 A = 2 J. According to a simple perturbation 
theory argument ^SiSL the polarization-averaged energy 
gap A = i J2 a Aq is, to a good approximation, the same 
as in the isotropic system with D = 0. The gap for exci- 
tations polarized along the x, y and z axes can then be 
written as: 



A, = A + 2D, 
A. = A,, = A — D. 



(2) 



The observable splitting D is proportional to the mi- 
croscopic anisotropy parameter D in the Hamiltonian. 
Numerical simulations indicate that for a uniform chain 
D « |.Di2L2& For isolated dimers one simply has D = D. 
For NTENP the gap energies can be estimated from high- 
field magnetization measurements of Narumi et al*^ The 
critical field at which the gap for one of the spin polar- 
izations is driven to zero by the Zeeman effect is given by 
g a [i B H { c a) = ^A^K;&&M For NTENP H { c x) = 9.3 T, 



H, 



( y) w H^ x) = 12.4 T, and g = 2.14. This gives 
A y « A z = 1.15 meV and A^ = 2.06 mcV. From this 
one gets A = 1.45 meV (A/ J w 0.43) and D « 0.3 mcV 

CD/J«0.1)4& _ 

The important energy scales A, hto^, v and D can be 
straightforwardly measured in inelastic neutron scatter- 
ing experiments, by mapping out the dispersion relation 



of magnetic excitations. They do not, however, carry 
any direct information on the level of dimerization in the 
system. 



C. Exploiting the lst-moment sum rule 

Additional insight can be drawn from an analysis of 
neutron scattering intensities of magnetic excitations. In 
fact, these intensities directly relate to the strengths of 
individual magnetic bonds. One way of extracting this 
information is by making use of the Hohenberg-Brinkman 
first moment sum rule for the magnetic dynamic struc- 
ture factor li^ 3 ^ This method has been previously suc- 
cessfully applied to the analysis of inelastic neutron scat- 
tering data on several occasions: for recent examples see 
Refs. I33U3 ll For the Hamiltonian Q the sum rule, an 
exact expression, can be written asiS 
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- ^2Ji sin 2 (qd 1 /2)(l - 8 a p){S ( 2 f S ( 2 f +1 ) 



- 2 J 2 sin 2 (qd 2 /2) (1 - 8 a0 )(S ( 2 f S ( 2 f_ l ) 



-2D{1 8 xf3 ) [2<(5f )) 2 ) + {{Sff) S(S + 1)] .(3) 

Here di and d 2 are real-space vectors chosen along the 
short and long bonds in the chains, respectively, and a 
and P label the coordinate axes: x, y and z. In practice it 
may be quite difficult to separately measure all three di- 
agonal components of S aa . Fortunately, for NTENP the 
ratio D/J is only about 10% (see discussion above), and, 
to a good approximation, the correlators (S^ S^j+i) anc ^ 
{S^j S^j-i) are independent of the subscript (3. The last 



term in Eq. scales as (D/J) 2 fRef. I26T and can be en- 
tirely neglected in our case. Under these assumptions the 
sum rule becomes: 



(huj)S aa {ci,uj)d{huj) 



-^ 1 sin 2 (qd 1 /2) - ^ 2 sin 2 (qd 2 /2). (4) 

The quantities E\ = ^(82^82^+1) and E 2 = 
J2(S2jS2j-i) are ground state exchange energies associ- 
ated with the strong and weak bonds, respectively. Due 
to the translational invariance, they do not depend on 
the site index j. 
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Eq. 0] directly relates the intensities measured in an in- 
elastic neutron scattering experiment to the modulation 
of exchange energy in the spin chains 



/(r. I. u.) 



s 



Ei — E2 
Ei + E2 



(5) 



While certainly not equivalent to 8, 5 is a very natu- 
ral measure of the magnitude of "dimerization" of the 
ground state. In the isotropic (Heisenberg) case of D = 
it can be directly expressed through S and the ground 
state energy E(S) = (H): 



2(1 + 5)S 2j S 2j+1 - 2(1 ~ 5)S 2j . 1 S 2j )/E(8) 



(6) 



For practical applications Eq. 0] can be further sim- 
plified if one assume the single mode approximation 
(SMA)^ 



5 M (q, W ) K 5 aa (q)i( W -^ a ). 



(7) 



The sum rule for a bond-alternating chain is then written 
as: 

<S QQ (q) ~~^-[Ei sin 2 (qd 1 /2) + £; 2 sin 2 (qd 2 /2)] . 

(8) 

For isolated dimers expression [7| is exact. For a uniform 
spin chain the SMA works remarkably well in most of the 
Brillouin zone, especially in the vicinity of the ID AF 
zone-center where the Haldane gap is observed ^^21 
Near the quantum critical point, where A vanishes, the 
SMA will fail entirely. However, for NTENP A/ J is sim- 
ilar to that in a uniform chain, and the SMA should still 
be reliable near the ID zone-center qo, qo(di + d 2 ) = 2tt. 
In this range the SMA dispersion relation can be written 
in the standard "relativistic" form: 
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'' sin 2 (qa) 



(9) 



FIG. 3: Contrast ration C(q) (Eq. for distinguishing 
ground state exchange energies based on inelastic neutron in- 
tensities as a function of momentum transfer in the (h, 0, 1) 
reciprocal-space plane of NTENP. The circled areas are re- 
gions of reciprocal space where most of the inelastic data were 
measured. 



that Ei and E 2 can not be distinguished based on the 
data collected at the ID AF zone-centers, where C(q) 
vanishes. This is rather unfortunate, since it is at these 
wave vectors that the dispersion w q is a minimum, and 
excitation intensity is maximized due to the l/w q fac- 
tor in Eq. [S] A high contrast is achieved away from the 
ID AF zone-centers, and at large momentum transfers. 
However, under these conditions the magnetic scattering 
is weakened by the 1 /w q coefficient and the effect of ionic 
magnetic form factors. Away from the AF zone-centers 
the applicability of the SMA also becomes questionable. 
Finally, the phonon background becomes progressively 
important at large |g| and interferes with the measure- 
ments. In our experiments we have found that a reason- 
able compromise between intensity, contrast and noise 
level can be achieved on either side of the h = 3 ID zone- 
center. Most of the data described in Section liV CI below 
were collected in that region of reciprocal space, repre- 
sented in Fig. |3| by the circled area in the lower-right. 



D. Application to NTENP 

The main experimental difficulty in using Eq. |H1 to 
estimate the ratio <5 in NTENP is the fact that the 
bond vectors di = 0.521a + 0.0246b - 0.0424c and 
d 2 = 0.479a - 0.0246b + 0.0424c are quite close in this 
material. Fig. |3 shows a grayscale and contour plot of 
the effective "contrast" ratio 

c( x sin 2 (qd 1 /2)-sin 2 (qd 2 /2) 
W sin 2 (qd 1 /2) + sin 2 (qd 2 /2) 1 ' 

as a function of momentum transfer in the (h, 0, 1) 
reciprocal-space plane in NTENP. This ratio is a good 
measure of our sensitivity to 5. Immediately one can see 



III. EXPERIMENTAL PROCEDURES 

Translucent dark-purple plate-like single crystal sam- 
ples of 50% deuterated NTENP were grown in aqueous 
solution. Four such crystals were co-aligned in one "su- 
persample" with a total mass of 1.5 g and mosaic spread 
of 2.5°. Inelastic neutron scattering experiments were 
performed at two different facilities. Lower-energy ex- 
citations (up to 3 meV energy transfer) were investi- 
gated using the SPINS cold-neutron 3-axis spectrome- 
ter installed at the NIST Center for Neutron Research. 
Neutrons with a fixed final energy Ef — 3.7 meV were 
utilized with a BeO filter after the sample. A Pyrolitic 
Graphite PG(002) monochronmator was used in com- 



5 



bination with a flat (Setup I) or horizontally focusing 
(Setup II) PG analyzers. Additional beam collimation 
was provided by the neutron guide and (open) — 80' — 
80' — (open) collimators (Setup II employed a radial post- 
sample collimator). Thermal-neutron studies were per- 
formed using the HB-1 3-axis spectrometer installed at 
the High Flux Isotope Reactor, Oak Ridge National Lab- 
oratory. The data were collected with Ei = 13.5 meV 
fixed-incident energy neutrons (Setup III). PG(002) re- 
flections were employed in both monochromator and an- 
alyzer. A PG filter was installed in front of the sample to 
eliminate higher-order beam contamination. Beam colli- 
mation was 48' - 40' - 40' - 240'. 

The sample being only partially deuterated led to a 
substantial geometry-dependent attenuation of the neu- 
tron beam due to incoherent scattering from hydrogen 
nuclei. This effect is equivalent to neutron absorption, 
and can be fully compensated for using the technique de- 
scribed in Ref.l38l For every inelastic scan measured, one 
performs a separate elastic scan to determine the effec- 
tive neutron transmission corrections. For every point of 
the "transmission" scan the sample rotation and scatter- 
ing angles are set exactly as in the inelastic scan. Unless 
a Bragg condition is accidentally satisfied in the sample, 
the main contribution to scattering in the "transmission" 
scan is due to incoherent elastic processes in the sample. 
To a good approximation, the corresponding cross sec- 
tion is isotropic and independent of neutron energy. The 
intensity detected in the "transmission" scan is there- 
fore directly proportional to the neutron transmission in 
the sample. Normalizing the original inelastic scan by 
the measured transmission correction not only gets rid 
of absorption effects, but also compensates for the any 
geometric corrections that occur when a large asymmet- 
ric sample rotates in a finite-size neutron beam in the 
course of the scan. In various scans measured in the 
present work the effective transmission coefficient varied 
by roughly a factor of 2 in the course of each scan or 
between different scans. 



IV. EXPERIMENTAL RESULTS 
A. Gap energies: constant-q data 

The gap energies and anisotropy splitting of the excita- 
tion triplet were accurately measured using cold-neutron 
Setups I and II. Energy scans collected at ID AF zone- 
centers (1,0, 0) (Fig.Ht) and (3, 0, 0) (Fig.gp) with Setup 
I show a single sharp peak at about 1.2 meV energy trans- 
fer. The background for these scans was measured, point- 
by-point at q = (1.5, 0, 0) and (2.5, 0, 0), respectively. As 
shown, the scans are corrected for transmission effects. 
The data in Fig. 0] were collected with the scattering 
vector q directed perpendicular to the (b, c) crystallo- 
graphic plane, and therefore represent fluctuations of y 
and z spin components of the triplet. The resolution of 
the present experiment is insufficient to unambiguously 
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FIG. 4: Inelastic scans collected at the (1,0,0) and (3,0,0) 
ID antiferromagnetic zone-centers in NTENP (symbols). The 
data were taken using a cold neutron 3-axis spectrometer with 
a flat analyzer. The solid lines are fits to the data as described 
in the text. 
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FIG. 5: Inelastic scan at the q = (1,0, 1) ID antiferromag- 
netic zone-center measured using a horizontally-focusing an- 
alyzer. The heavy solid line is a fit to the data as described 
in the text. The light-grey and dark- grey shaded areas repre- 
sent partial contributions of gap excitations polarized parallel 
and perpendicular to the crystallographic (6, c) plane, respec- 
tively. 



resolve the gaps for y- and z-polarized modes. To detect 
the x-axis spin fluctuations we performed additional mea- 
surements with a large momentum transfer perpendicular 
to the chains. In the focusing-analyzer mode (Setup II) 
the scattering vector was at all times maintained on the 
(1,0, 1) reciprocal-space rod. For each energy transfer, 
the transverse momentum transfer / was chosen to have 
the sample chain axis parallel to the scattered neutron 



6 
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FIG. 6: Constant- E scans measured in NTENP in the vicin- 
ity of the for (1,0,0)AF zone-center (symbols). The dashed 
and dotted lines in (a) are profiles simulated assuming a cross 
section function as given by Eqs.[HJ[j3witli S = — 1 and 8 = 1, 
respectively. The sold line in (a) is a global fit to several 
scans, as described in the text. The solid lines in (b) and (c) 
are fits to single scans. 



beam. This geometry optimizes ^-resolution along the 
chains. The resulting scan is shown in Fig. EI Due to the 
intrinsic polarization-dependence of the magnetic scat- 
tering cross-section, the largest contribution to this scan 
is from y and cc-polarized excitations. The new feature 
observed at about 1.8 meV energy transfer was thus at- 
tributed to s-axis spin fluctuations. Note that the lower- 
energy peak in Fig. [S] appears at a slightly lower energy 
than in the scan taken with the scattering vector parallel 
to the chain axis. This behavior is likely due to a weak 
dispersion of excitations perpendicular to the chains. 

The constant-Q scans were analyzed assuming a sim- 
ple single-mode cross section as given by Eqs. |S] and 
|5J The spin wave velocity was fixed at v = 8.6 meV, 
as separately determined from the analysis of constant- 
E scans described below. The usual polarization factors 
for unpolarized neutrons determined the relative intensi- 
ties of the in-plane and out-of-plane spin fluctuations. 
In addition, a magnetic form factor for Ni 2+ was ex- 
plicitly included in the cross section function. The re- 
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FIG. 7: Measured dispersion relation of magnetic excitations 
in NTENP (symbols). Lines are the result of a global fit of 
a model cross section to the neutron data as described in the 
text. 



suiting model for the dynamic structure factor was nu- 
merically convoluted with the spectrometer resolution 
function, calculated in the Cooper-Nathans approxima- 
tion. The only adjustable parameters of the model were 
the gap energies A^, A y and A z and an overall inten- 
sity scaling factor. A very good fit to the scans with 
zero transverse momentum transfer was obtained with 
A y — A z — 1.16(0.01), as is shown in solid lines in 
Fig. 0] Slightly smaller energies of the doublet were ob- 
tained by analyzing the data shown in Fig. |3J Here the 
gap energies were found to be A y = A z — 1.07(0.01) 
and A x = 1.91(0.02) meV. These values correspond to 
A = 1.35(2) meV and D = 0.28(2) meV, in good agree- 
ment with bulk measurements of Narumi et alm^ 



B. Dispersion relation: constant-i? data 

The dispersion relation for single-mode excitations in 
NTENP was measured in constant-i? scans. In order 
to optimize magnetic intensities (form factor) these data 
were collected around the (1,0,0) AF zone-center (left 
circled area in Fig. |3J). Typical scans are plotted in open 
symbols in Fig. El To take into account resolution ef- 
fects, but to avoid constraining the dispersion relation 
to the postulated sinusoidal form, each scan was sepa- 
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FIG. 8: Constant--^ scans measured in NTENP using a 
thermal-neutron setup in the vicinity of the (3, 0, —2.5) AF 
zone-center (symbols). All lines are as in Fig. |SJa). 



rately analyzed using the model cross section described 
above. The adjustable parameters for each scan were an 
intensity prefactor and the spin wave velocity v. The gap 
energies were fixed at the values determined using cold 
neutrons (see above) . Typical fits are shown in Fig. EIb) 
and (c) in solid lines. The dispersion relation deduced 
from such fits to individual scans is plotted in open sym- 
bols in Fig. 



C. Exchange energy modulation and global fits 

All scans described above are fairly insensitive to the 
distribution of exchange energies in the chains, due to a 
small contrast ratio at the particular wave vectors. To 
illustrate this, in Fig. Eta) we have plotted peak profiles 
simulated using Eqs. [2HSI and assuming E\ = (dashed 
line) or E<z = (dotted line). To within the accuracy 
of our measurements, the two profiles are almost iden- 
tical. A much better contrast (although a significantly 
smaller intensity) was achieved in constant- E scans col- 
lected along the (h, 0, —2.5) reciprocal-space rods around 
h = 3 (right circled area in Fig. [2). These data are plot- 
ted in symbols in Fig. [H] A similar constant-^ scan was 
measured using cold neutrons and Setup I, and is plotted 
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FIG. 9: Constant-iJ scan measured in NTENP using a cold- 
neutron setup in the vicinity of the (3, 0, —2.5) AF zone-center 
(symbols). All lines are as in Fig.|SJa). 



if Fig. El In both figures the dashed and dotted lines are 
simulations for E 2 — or E\ — 0, respectively, assuming 
v = 8.6 meV. From these data the strong dimerization in 
NTENP becomes apparent: the ground state exchange 
energy is primarily concentrated on the shorter bonds. 

To quantify this observation we performed a global fit 
to the data collected in constant-i? scans measured for 
energy transfers up to 5 meV with setups I and III. The 
adjustable parameters were two intensity scaling factors 
(one for each setup), the spin wave velocity v and the 
bond-alternation parameter S. Very good fits are ob- 
tained with v = 8.6(1) meV and 5 = 0.42(+0.2, -0.1). 
The large asymmetric error bar on 5 is unavoidable due 
to the technical difficulties associated with low intensity, 
small contrast ratio and transmission corrections. 

The simple sinusoidal dispersion curves postulated in 
Eq. El is plotted in lines in Fig. using the experimen- 
tally determined gap energies and spin wave velocity. At 
higher energies the experimental data points clearly lie 
below these curves. While it is difficult to extrapolate 
the experimental result to the zone-boundary, one can 
roughly estimate Hlozb to be smaller than v by about 
1 meV. 

All physical parameters obtained for NTENP in our 
neutron scattering studies are summarized in Table [I]m 
comparison with those obtained by Narumi et al. using 
bulk techniques 14 and to known exact results for uniform 
and fully dimerized isotropic 5 = 1 chains. 



V. NUMERICAL CALCULATIONS 

In order to relate the measured exchange energy mod- 
ulation parameter S to the alternation of exchange con- 
stants S, we performed a numerical study of the model 
Hamiltonian l(T]). assuming a vanishing anisotropy D = 0. 
The ground state energy E{5) was computed as a func- 
tion of 5 using the density matrix renormalization group 
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TABLE I: Physical parameters for the alternating S — 1 chains in NTENP in comparison with those for uniform and fully 
dimerized isotropic S — 1 spin chains. 

NTENP" NTENP 6 Uniform chain Isolated dimers 



J 3.4 meV 

A 1.45 meV w 0.42J 1.35(2) meV « 0.40(1) J 0.41J 2J 

D 0.3 meV » 0.1 J 0.28(2) meV » 0.083(6) J 2^27,28 D 

i! 8.6(1) meV» 2.5(1)J 2.5^ 

fitJzB » 7.5(5) meV= 2.2(1) J 2J 

S = E i- E 2 _ o 42 1 

El +E 9 , ^ ^ 

5 = i^T 1 °- 40 or ~ C 0.37(1) or 0.06(2) d 1 



0.30(0.05) e 
0.24(-0.04,+0.08) / 



"Bulk measurements: Ref. ll4lisl 

'Neutron scattering: this work, assuming J = 3.4 meV. 

c From A/ J with A deduced from high-field magnetization data 
and numerical results of Ref . 0. 

^From A/ J with directly measured A and numerical results of 
Ref.0 

e From Hujzb/J with directly measured Hlozb and numerical re- 
sults of Ref. 2a 

■'From directly measured 5 and numerical results of this work. 



(DMRG) method 39 i 40 i 41 for a chain of 32 spins with pe- 
riodic boundary conditions. The parameter 5 was then 
obtained using Eq-El These results are plotted in Fig. 1101 
The computed curve is monotonic and quite smooth in 
the studied domain of parameter space. Crossing the 
quantum-critical point at S c w 0.26 is marked only by 
a weak point-of-inflection-type anomaly. We recall that 
at S c the correlation length diverges. However, for our 
particular calculation, this does not appear to be a prob- 
lem. The observable ((HJ is very smooth and is not subject 
to significant finite-size corrections. It was verified that 
doubling of the system size changed the results shown in 
Fig-EUby less than one percent. 

To justify the use of the isotropic model for computing 
5 as a function of S, we also studied the behavior of the 
expectation value of the anisotropy term in the Hamil- 
tonian^ It was found that that its contribution to the 
ground state energy is practically independent of 5 and 
is equal to about 2 /3-D per spin. 



VI. DISCUSSION 




cHVJ 2 )/(j 1+ j 2 ) 

FIG. 10: Alternation of exchange energies 8 — (Ei — 
E2)/{Ei + E2) as a function of the alternation of exchange 
constants 8, computed for a 32-site bond-alternating S = 1 
chain with periodic boundary conditions. The quantum phase 
transition at S c corresponds to an inflection point on the cal- 
culated curve. 



Our neutron results present new opportunities to de- 
termine the magnitude of bond-alternation in NTENP 
and to place this material on the phase diagram. First of 
all, S can be estimated from the gap energies. Using the 
numerical results of Ref. and assuming J =3.4 meV 14 , 
for A = 1.35(2) meV one gets S = 0.37(1), assuming 
NTENP is on the dimerized side of the phase diagram. 
The alternative is S = 0.06(2) in the Haldane phase, and 
the distinction can not be made based on gap measure- 



ments alone. 

The numerical results of the previous section allow 
us to independently estimate S based the measured ex- 
change energy alternation. The experimental value 6 = 
0.42(+0.2, -0.1) corresponds to5w 0.24(-0.04, +0.08). 
This value is within 1.5tr of the estimate based on the 
measured gap energy, with the assumption that NTENP 
is in the dimerized phase. Assuming that NTENP was in 
the Haldane phase with S = 0.06 would imply i5 ~ 0.05, in 
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considerably worse agreement with experiment. The er- 
ror bars associated with our measurements of S are rather 
large, and the method itself relies on the uncontrolled 
single-mode approximation. However, when combined 
with gap measurements, these data cleraly confirm that 
NTENP has a dimerized, rather than Haldane ground 
state. The same conclusion can be reached by compar- 
ing the measured zone-boundary energy to numerical re- 
sults of Ref. |2^, according to which our experimental es- 
timate Tilozb/ J ~ 2.2(1) corresponds, unambiguously, to 

5 = 0.3(0.05). 

It is important to emphasize that our combination 
experimental approach allows to uniquely determine 6 
and decide whether the ground state is dimerized or 
not, based on measured properties of the undoped ma- 
terial. This eliminates a possible ambiguity associated 
with guessing the ground state from the behavior of 
non-magnetic impurities, as was done in Ref. 14. The 
problem is that the ground states of chain- fragments in 
the impurity-doped system are different from the ground 
state of defect-free chains. In particular, for the Haldane 
phase (5 < S c ) there are effective interactions between 
the two liberated S = 1/2 spins on the ends of each 
chain fragment 42iM44 Depending on the parity of the 
number of magnetic sites in the fragment, the interac- 
tion is antiferro- or ferromagnetic^ In the former case, 
the fragment has a non-magnetic 5 = ground state 
and does not produce any ESR signal. In odd-length 
fragments, however, end-chain spins combine to form an 
5=1 triplet. At sufficiently low temperature any ESR 
experiment will observe S = 1 (rather than 5 = 1/2) 
degrees of freedom, just like in the dimerized phase, for 

6 > S c . 



VII. CONCLUSION 

In summary, we have performed a combination of in- 
elastic neutron scattering measurements to determine the 
key characteristics of the bond-alternating 5 = 1 quan- 
tum spin chain compound NTENP. The results unam- 
biguously indicate a strong but incomplete dimerization 
of the ground state in this material. Future studies, per- 
haps using larger and fully deuterated single crystal sam- 
ples, will concentrate on features of the spectrum that are 
beyond the simplified single-mode picture. 
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